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Introduction 



This note describes the results of an exploratory data- 
analytic study of sea surface temperatures and wind speeds meas- 
ured on Oceanic Weather Ship PAPA. The data consist of 15 years 
(1955-1969) of measurements taken every three hours. Two of the 
measurements recorded were of wind speed and sea surface tempera- 
ture. other measurements taken include air temperature and wind 
direction . 

In Section 2 the results of a spectral analysis of the 15 
years of sea surtace temperatures are described. It indicates 
that sea surface temperatures exhibit noticeaole 1 year, b month, 
24 hour, and 12 hour cyclic components. The 24 hour and 12 hour 
components are modulated by yearly and seasonal effects. 

In Section 3 the results of an analysis of average daily sea 
surface temperatures and average daily wind speeds during the 
Spring and Summer for three years 1964-66 are reported. After the 
two series were detrended, it was found that positive changes 
between days t and t-1 in the inverse residual sea surface 
temperature are associated with large residual wind speeds on day 
t and large residual sea surface temperatures on day t-1 . 

This association can be explained by the behavior of the oceanic 
mixed layer depth. 
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1 . Sea Surface Temperatures 



In this section we describe a spectral analysis of sea sur- 
face temperatures taken on Ship PAPA. The data consist of 15 
years (1955-1969) of measurements taken every three hours. The 
5 leap days were removed leaving a data set of length 43,800 
readings . 

Figure 1 shows the data set. The most apparent feature of 
sea surface temperatures is the strong yearly cycle. Another 
feature is the change in the measurement reporting procedure 
which occurs 8 years into the record. The first 8 years of data 
appear to be more discrete than the last 7 years. Prior to 1963 
the data were recorded with a 1° Fahrenheit resolution and then 
converted to degrees Celsius. Starting in 1963 the data were 
recorded with a 0.1° C resolution. The data record also appears 
to contain erroneous measurements and interpolations. 

The 2n (normalized periodogram) of the entire record of sea 
surface temperatures appears in Figures 2-5; (cf. Cox and Lewis 
[1966] p. 99). Figure 2 shows the entire Jin periodogram; the 
horizontal line is at the Jin 95% quantile of the maximum of 
21,900 independent unit exponentials. Figure 3 shows the first 
hundred values of the Jin (periodogram). The two peaks occur at 
p = 15 and 30 which correspond to cycles of 1 year and 6 months 
respectively. Figure 4 shows the values of the Jin (periodogram) 
for p-values around the peak at p = 5475 which is the 24 hour 
cycle. Notice the local peaks at p = 5490 and 5460 which are 
1 year side lobes. A second pair of side lobes occur at p = 5408 
and 5542 which correspond to 81.7 day side lobes or seasonal side 
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lobes. Figure 5 gives the £n ( per iodogram ) for values of p 
around the peak at 1U95U, the 12 hour cycle. Again there are 
side lobes at p = 1U945, 1U965 (1 year) and p = 10882, 11018 
(80.5 days or seasonal). 

The side lobes indicate that sea surface temperatures are 
modulated by yearly and seasonal effects. It suggests that the 
sea surface temperatures in the same season in different years are 
statistically different; (ct. Bloomfield [19 7b] pp . 99-100). 

R. w. Garwood has provided the following explanation tor the 
year to year variability. The variability in sea surface tem- 
peratures in the same season of different years is due to the 
variability of the mixed layer depth in the ocean. The mixed 
layer depth is a function of surface heating of the ocean and 
ocean mixing by turbulence due to storms. The times between 
passage of, and the strength of, the storms varies from year to 
year. The surface heating is a function of the season of the 
year . 

In summary, the In ( penodogram ) of 15 years of sea surface 
temperatures indicate that the record has noticeable 1 year, 
b month, 24 hour, and 12 hour cyclic components. The 24 hour and 
12 hour cycle components are modulated by year and seasonal 
ef tects . 
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2 . Associations between average daily Sea burface Temperatures 
and Average Pail/ wind Speed 

In this section we will report on a pilot study of associa- 
tions between averaye daily sea surface temperature and averaye 
daily wind speeds duriny the Spring and Summer (Julian days 
91-272) tor the years 1964 (year 10), 196b (year 11), and 1966 

(year 12 ) . 

P. A. W. Lewis and his associates have done a spectral anal- 
ysis of 15 years (1955-1969) of Ship P wind speed data which 
was measured every three hours ( Lewis 1 1963 ]) . The analysis sug- 
gests that the wind speed data has cyclic components of 1 year, 

6 month, 12 hours and 6 hours. 

Associations between the average daily sea surface tempera- 
tures and average daily wind speeds may be confounded by the non- 
stationarity of the two series. Thus, the two series were de- 
trended. In order to detrend the series in the same manner, the 
average daily sea surface temperatures and average daily wind 
speeds for Julian days 91-272 for each of years 1U-12 were put 
into two-way tables having 26 rows (weeks) and 7 columns (days). 
The two-way tables were median polished (cf. McNeil [1977J) to 
obtain expressions of the form Data = typical value + day effect 
+ week effect + residual. The residuals from the median polish 
were taken as the detrended series. Since both averaye daily sea 
surface temperatures and averaye daily wind speeds were treated 
in the same manner by this procedure, the residual series are 
comparable. figure 6 (respectively 7) shows a plot of the median 
polish residuals of averaye daily sea surface temperatures 
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(respectively wind speeds) for year 12. The plots for the other 
years are similar. The plots exhibit no apparent nonstat ionarity 
and show the typical pattern of many small values due to the use 
of median polish. The periodogram (not shown) of the residual 
sea surface temperature for year 12 shows a significant peak at 7 
days which may be due to the fact that the median polish was done 
on a two-way table with rows of length 7 days. The periodograms 
of residual sea surface temperatures for the other two years do 
not show this effect. 

The first 50 serial correlations were computed for each 
residual series. There is some correlation in the residual 
series particularly in the residual average sea surface tempera- 
tures in year 12. There appears to be a pattern in the correla- 
tions again suggesting possible nonstat ionarity introduced by the 
median polish. However, the correlation was judged not to be 
substantial enough to affect the results of this exploratory 
study . 

Let W ( t ) (respectively S(t)) denote the residual from 
median polish of the average daily wind speed (respectively sea 
surface temperature) on day t . Put 



H( t ) 



1 

S(t)+1 ' 



t = 91, ... ,272. 



where the constant 1 in the denominator was arbitrarily chosen 
to ensure that H(t) is finite. The surface mixed layer depth 
of the ocean during the Spring and Summer is related to the in- 
verse of sea surface temperature; (Garwood [1983]). Thus, H(t) 
should be related to mixed layer depth. 
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Let 



D ( t ) = H(t) - H(t-l) , t = 92, ..., 272 

and 

X ( t ) = £n (S ( t ) + 1) t = 91, ..., 272. 

D(t) is related to changes in mixed layer depth. We will con- 
centrate on associations between D(t) and W(t) and X(t-l). 

For each year the quantiles of residual sea surface tempera- 
ture fb(t) ; t = 91, ...,271} were used to categorize 

{ ( D ( t ) , W(t)), t = 92, ...,272} in the following manner. Category I 
contains those (D(t), W(t)) such that S(t-l) is less than or 
equal to the lower quartile of {S(t); t = 9 1 , . . . , 272 } ; Category II 
contains those data such that S(t-l) is greater than the lower 
quartile but less than or equal to the median; Category III contains 
those data such that S(t-l) is greater than the median but less 
than the upper quartile; and Category IV contains the remainder 
of the data. In brief. Category I is "very low"; Category II is 
"low"; Category III is "high"; and Category IV is "very high" as 
ordered by sea surface temperature on the previous day. 

Figures d-lu show plots of residual wind speed W(t) versus 
D(t) by residual sea surface temperature category for each of the 
three years. For each year the plot in the upper (lower) left hand 
corner corresponds to Category I (Category III); the plot in the 
upper (lower) right hand corner corresponds to Category II 
(Category IV) . 

All of the plots suggest that larger residual wind speeds are 
associated with larger values in the differences D(t). Further the 
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is no 



association looks roughly linear. There 
that the association changes for differen 
categories . 

For each year the quartiles of {w(t 
used to categorize {(D(t), X(t-l)), t = 
manner. The categorical plots of X(t-l) 
in Figures 11-13, the arrangement of the 
same as for Figures 8-10. Once again the 
larger values of D(t) are associated wi 
X(t-l). The association appears roughly 1 
strong indication that the association is 
categories . 

Regressions of the form 

D(t ) = a + b W(t ) + cX(t-l) 



t sea 



strong indication 
surface temperature 



) ; t = 91,. . . ,272} were 
92, ... 27 2} in the same 
versus D(t) are shown 
plots by category is the 
plots indicate that 
th larger values of 
inear. Again there is no 
different in different 



t = 92, . . . , 27 2 



were fit to the data in each of the three yea 
cedures. One is least-squares, a classical p 
sensitive to possible outlying values in the 
and Tukey [1977]). The second is a biweight 
Mosteller and Tukey [1977] pp . 205) which is 
outlying values. The estimated coefficients 
Table 1. Plots of the residuals from the reg 
dieted D ( t ) showed little structure. This 
regression has removed most that is systemati 
the data. 

Regressions containing an interaction te 
D ( t ) = a + b W(t) + cX(t-l) + d(X(t 
were also fit. The estimated values appear i 



rs using two pro- 
rocedure which is 
data; (cf. Mostelle 
procedure (cf. 
less sensitive to 
can be found in 
ression versus pre- 
suggests that the 
c or explainable in 

rm of the form 

-1) x w(t) ) 

n Table 2. 
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Table 1 



Estimated Values for the Coefficients of 
D ( t ) = a + b W ( t ) + cX( t-1 ) 



Year 


Least 


Squares Estimates 
(Std. Error) 


Biweiyht Estimates 




a 


b 


c 


a 


b 


c 


10 


0.0141 
( .017) 


0 . 202 7 
( .048 ) 


0.6510 
( .090 ) 


-0.0081 


0.1267 


0.8510 


11 


0.0144 

(0.027) 


0.1325 

(0.074) 


0.8969 
( 0.108 ) 


-0.0214 


0.1319 


0.6120 


12 


0.0125 

(0.024) 


0.1475 

(0.063) 


0.7368 

(0.092) 


-0.0179 


0.1368 


0.4254 



Table 2 

Estimated Values for the Coefficients of 
D(t) = a + b W(t) + cX(t-l) + d(X(t-l) x w(t)) 



Year 


Least squares Estimates 
(Std. Error) 


Biweiyht Estimates 






a 


b 


c 


d 


a 


b 


c 


d 


10 


0.0188 
( 0.017 ) 


0 . 2 1 7 5 
(0.048 ) 


0.6298 
( 0 .090 ) 


-0 .6047 
( 0.311 ) 


0.0104 


0.1144 


0.5202 


0.2b 


11 


0.0153 
(0.027 ) 


0.1504 
( 0.0 7 5 ) 


0.8328 

(0.115) 


0.4474 
( 0.272) 


-0.0210 


0.1346 


0.6114 


-0.09 


12 


0.0314 
( 0.024 ) 


0.1933 

(0.062) 


0.7193 
(0.089 ) 


-0 .9222 
(0.236) 


-0.0202 


0.1312 


0.4023 


0 . 20 
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The regressions indicate that larger differences in the 
inverse residual sea surface temperature, D(t) , are associated 
with larger residual wind speeds, W(t) and larger £n (residual sea 
surface temperature) the day before, S(t-l) . There appears to 
be no significant interaction term of the form W(t) x X(t-l) . 

The standard errors for the least squares estimates should be 
viewed with caution because of the correlation in the series 
(S(t)} and {W(t )} . There is no suggestion that the association is 

different in the different years. 

An explanation for the relationship suggested by R. W. Garwood 
is as follows. The mixed layer depth of the ocean during the Spring 
and Summer is related to the inverse of sea surface temperature. 

Thus D ( t ) should be related to the change in mixed layer depth of 
the ocean between days t and t-1 . Further 

X(t-l) = £n(l + S(t-l) ) is related to the mixed layer depth on 
day t-1 . Thus the above association suggests that deepening 
(shallowing) of the ocean mixed layer depth during the Spring and 
Summer (Julian days 92-272) is associated with larger (smaller) wind 
speeds and shallower (deeper) mixed layer depths the day before. 

During the Spring and Summer the mixed layer depth of the 
ocean is a function of solar warming and turbulant mixing attri- 
butable to storms. During periods of light (strong) winds the 
mixed layer becomes shallower (deeper). However, the effect of a 
storm on the mixed layer depth depends on "how deep the mixed layer 
is already". A storm of the same strength will cause a larger drop 
in sea surface temperature if the surface mixed layer is shallow 
than if it is deep. 
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The association found here agrees with a simpler related 
association found by Elsberry and Raney (ly78) also using bhip P 



data. They defined sus 
wind forcing as events 
change for each event, 
increase (decrease) in 
average wind speeds. 



tained periods of above or 
and calculated the surface 
They found an association 
sea surface temperature and 



below normal 
temperature 
between an 
below (above) 
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Conclus ions 



i . 

(1) Spectral Analyses of 15 year records of sea surface 
temperatures ana wind speeds from Ship P showed strony yearly, 
seasonal and time of day effects. 

The Spring-Summer periods of three years 1964-1966 were chosen for 
further study. 

(2) After average daily sea surface temperatures and wind 
speeds for the period were detrended and the residuals examined, 
it was found that large changes between days t and t-1 in the 
inverse residual sea surface temperature are associated with 
large wind speeds on day t and large sea surface temperatures 
on day t-1 . The association can be explained by the behavior 
of the oceanic mixed layer depth. 
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